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The current 3 ct discrepancy between experiment and Standard Model predictions for 7r° —>■ e''"e“ 
is reconsidered nsing the Fade Theory for bivariate functions, the Canterbury approximants. This 
method provides a model-independent data-driven approximation to the decay as soon as experimen¬ 
tal data for the donbly virtnal 7r° transition form factor are available. It also implements the correct 
QCD constraints of the form factor both at low- and high-energies. We reassess the Standard Model 
result including, for the first time, a systematic error. Our result, BR(7r° —>■ e^e“) = 6.23(5) x 10“®, 
still represents a discrepancy larger than 2cr, unsurmountable with our present knowledge of the 
Standard Model, and would claim New Physics if the experimental result is confirmed by a new 
measurement. Our method also provides the adequate tool to extract the doubly virtual form factor 
from experimental data in a straightforward manner. This measurement would further shrink our 
error and establish once and for all the New Physics nature of the discrepancy. In addition, we 
remark the challenge this discrepancy poses in the evaluation of the hadronic light-by-light scatter¬ 
ing contribution to the {g — 2)^, specially confronted with the foreseen accuracy of the forthcoming 
(g — 2)fj, experiments. 
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I. INTRODUCTION 


Pseudoscalar decays into lepton pairs provide a unique 
environment for testing our knowledge of QCD. As such 
decays are driven by a loop process, they encode, at once, 
low and high energies. For the 7r° decay, the process 
(neglecting electroweak corrections) proceeds (Fig. 
through the 7r° —)■ 7*7* anomalous vertex [T], with the 
photons linked by a lepton line. The loop does not di¬ 
verge due to the presence of the tt^ transition form factor 
(TFF) on the anomalous vertex, the , (q — k)'^) 

with {q — k)"^ the photon virtualities. The TFF 



FIG. 1: Feynman Diagram for tt® —>■ e^e process. 

with —q^= cannot be calculated from first princi¬ 
ples, and only the = 0 and > 00 limits are known 
in terms of the axial anomaly in the chiral limit of QCD 
[T] and perturbative QCD [ 5 ], respectively. It has been 
customary to introduce a model for the TFF to interpo¬ 
late between both regimes through the use of vector me¬ 
son dominance ideas [SE] with different assumptions for 
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describing the doubly virtual behavior [ 5 ]. The discrep¬ 
ancy among different choices there reflects the model- 
dependency of that procedure (see, for example, the dis¬ 
cussion in Ref. [ 5 ]). In order to avoid this undesired 
model dependency, two main strategies have been devel¬ 
oped so far based on the observed low-energy nature of 
the process (c.f. Eq. ^ and discussion afterwards). The 
first one involved the use of chiral perturbation theory 
(xPT) to calculate the diagram in Fig. (see [6] and ref¬ 
erences therein), where the low-energy constants of the 
theory would encode the high-energy effects from the pro¬ 
cess which are common, to first approximation, to tt^ and 
T] allowing to predict the 7r° decay in terms of the mea¬ 
sured g one [ 7 ]. The second one involved the use 

of a phenomenological parametrization of the TFF fitted 
to experimental data [ 5 ]. The lack of double virtual data 
for the double TFF however has been complemented so 
far with additional high-energy QCD constraints, which 
might distort the low-energy region fundamental in this 
process. 


In this work, after briefly discussing the 


+0— 


e^e 


process in Sec. |TTj we explore for the first time the role 
of the bivariate Fade approximants, the so-called Can¬ 
terbury approximants (CA) [ 5 ] in Sec. Ill which will 
allow us to describe from the low-energies and in a 
model-independent data-driven approach the doubly vir¬ 
tual TFF driving this rare decay. CAs will allow us to 
evaluate the impact of the QCD high-energy tail in a 
consistent manner as well. We remark that the CA’s 
technique can be generalized straightforwardly to other 
processes involving analytic Stieltjes and meromorphic 
functions with two variables. In the present work we 
exemplify its usage in a theoretically challenging pro¬ 
cess. Later on, in Sec. |IV| we will consider the impact 
of the KTeV measurement on the HLBL contribution to 






{g — 2)e_^ using the CAs. A potential new physics sce¬ 
narios are considered in Sec. El and the main results of 
our work are collected in the Conclusions. 
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where /3e = (1 — is the outgoing lepton 

velocity and A(m^o) is given by the loop integral 


II. THE 7 r° ^ e+e" DECAY 


The most accurate measurement of the tt® —>■ e“'"e“ 
was performed by the KTeV Collaboration at Fermi- 
lab through the observation of almost 800 7 r° —>■ e+e” 
events [9] and yielded BR( 7 r° —>• e+e“) = (7.48 ± 0.29 ± 
0.25) X 10“® after removing the final state radiative cor¬ 
rections (RC) [TU], where the first error referred to statis¬ 
tics and the second to the total systematics. 

The normalized BR( 7 r° —?> e^e“) is defined as 


BR( 7 r° -» e+e ) QemWe 

BR(7r‘^ —>■ 77) \ Trm^o 
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( 1 ) 


A{q^) = 2i 


d'^fc , {q - kf) 

7^2 q^k^[q — kY[(jp — k)"^ — ml) 


and encodes all the hadronic effects through the normal¬ 
ized TFF ,{q — k)"^) (i.e. ^^ 077 ( 0 , 0 ) = 1). 

Even without any information about the TFF, Cutcosky 
rules may be used to extract its imaginary part, which 
provides the well-known unitary bound discussed by 
Drell [TT], BR(7rO ^ e+e") > BR“““^’'y(7rO ^ e+e") = 
4.69 X 10“®, which is a model-independent result. 

The presence of the photon propagators (cf. Fig. 
implies the kernel of the loop integral ([^ to be peaked 
at very low energies of around the electron mass as is 
shown in Fig. The kernel can be expanded in terms of 
me/m^o as well as mg/A and m^o/A, being A the cut-off 
of the loop integral, or the hadronic scale driven by the 
TFF. Then, Eq. ([^ reads ISE]: 


A(m^o 


ITT 1 

^ We ^ We 




Pe-l 
1 + /3e 


where L = In and terms of O In 

as well as C>(^^,^ln^^, and C>(^y,yin^j 
have been neglected. 

The integral in Eq. (§, Fig.§ produces a negative 
contribution, diminishing the result. Omitting such con¬ 
tribution, Eq. (|^ would result in 19 x 10“® for the BR. 

Recently, the authors of resummed the power cor¬ 
rections using the Mellin-Barnes technique and found 
them of 0(1%) [H]. Then, using a Vector Meson Dom¬ 
inance for the TFF, they found BR( 7 r° —>■ e+e”) = 
6.2(1) X 10“® OEU, 3.2ct off the KTeV result. 

Such discrepancy demands further explanations pro¬ 
vided that future experiments (for example, the NA48/2 
or NA62 kaon experiments at the CERN SPS which have 
demonstrated the possibility to perform precision tt® de¬ 
cay physics |18j l would confirm the current measurement. 
Three research lines can be conceived: a reevaluation of 
the radiative corrections, an improved parameterization 
of the doubly virtual TFF, or a new mechanism within 
physics beyond the SM [HIITS] . 

In Ref. |T7j , the radiative corrections used by the KTeV 
based on Bergstrom’s work m were reconsidered. At 
that time, Bergstrom considered the two-loop QED ra¬ 
diative correction to the decay in the soft-photon approx¬ 
imation together with an inclusion of a certain cut-off for 
the loop diagrams. He also considered the role of the 
Dalitz decay and its interference as a source of experi¬ 
mental error. The authors of noticed that m ne¬ 
glected a class of subleading diagrams, which due to par¬ 
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Q‘ 




(3) 


ticular cancelations among the dominant ones, turned 
out to be dominant. Later on, they also studied the 
role of the soft-photon approximation finding it accurate 
enough [TB]. Ref. m suggests that the radiative correc¬ 
tions represented a —13% effect, so increasing the value 
measured by KTeV (see [9] for details). The reanaly¬ 
sis of Refs. EZlllH] suggested, however, that including 
the subleading diagrams the RC would decrease down to 
— 6 %, implying a smaller BR after RC are removed. With 
such considerations, the new KTeV value would result in 
BR"KTeV” (tt® —t e+e“) = (6.87 ± 0.36) x 10“®, closer to 
the SM value at about 2cr. 


As we stated before, we investigate here the role of the 
TFF on such decay. Given the accuracy of Eq. ([^, it is 
safe to conclude that the main contribution to the loop 
integral ([^ happens at low space-like energies (Fig. |^. 
Therefore, a precise description of the doubly virtual 
TFF at low space-like energies is an essential starting- 
point for an accurate prediction, an observation over¬ 
looked so far. Moreover, we noticed that the factor¬ 
ization approximation for the TFF, i.e., F{Ql,Ql) = 
F{Ql,0) X F(0,Q2)) can induce large effects. For our 
study, we consider the reconstruction of the TFF of dou¬ 
ble virtuality by a data driven approach. That method, 
model independent, is based on the theory of Fade ap- 
proximants (PA) [TB] extended to the double virtual case, 
the Canterbury approximants [5]. 
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FIG. 2: Kernel of the integral in Eq. for different energy ranges. 


Each band labels the partial contribution to the integral. 


III. CANTERBURY APPROXIMANTS 

Canterbury approximants (CA) are rational approxi- 
mants defined from double power series in variables Qi 
and Q 2 with the following properties [ 8 ]: 

1 . CA are symmetric with respect to Qi and Qi- 

2. CA reduce to the PA if QfiQ^) = 0. 


CA, defined as C^{x,y) = 0; extend PAs 

from one to two variables. The coefficients of Rn{x, y) — 
ihj) € and QN{x,y) = Yn^hyx'^y^, 
{i,j) € "D are determined by the accuracy-through-order 
conditions up to order 0{x ^^and 
with a = 0 , 1 ,... 2 A, i.e., given f{x, y) = Yij Cijx^y^, 

a /3 

'y ^ y ^ — (tt, /?) € Af 

i=0 j=0 


3. CA enjoy homographic invariance under argument 
transformations. 

4. CA are unique. 


and ( 60,0 = 1 as part of the definition) 

min{a,N) mm{0,N) 

E E h,jCa-i,p-j = Ofor (a,/3) e £, (a,/3) ^ N 

i—0 j—0 


5. CA satisfy the low-energy constraints by definition. 

6 . CA easily accommodate the high-energy QCD con¬ 
straints. 

Even though generalizations of CA are known [20112^] , 
they will reduce to CA for the 7 r° e'^e~, represent¬ 

ing the simplest and most efficient model-independent 
approach. On top, the convergence to symmetric two- 
variable Stieltjes functions is known to exist [8l|24]. In 
our case of study, the TFF, the first property follows triv¬ 
ially from Bose symmetry. On top, when one of the pho¬ 
ton virtualities is zero, the resulting TFF can be safely 
approximated by a standard PA with great success [25l - 
HZ], satisfying the second CA property. The third prop¬ 
erty is used in [H] together with the accuracy-through- 
order conditions to explore the convergence of the CA. 
As it is shown below, properties 4), 5) and 6 ) emerge 
naturally by construction. The important remark here 
is that the TFF of doubly virtuality PSIET] is a bivari¬ 
ate Stieltjes function [21] for which the CA’s convergence 
theorem from [El ED apply. 


where dim(£’) = dim(A/') + dim(I?) — 1 . With x = Qi 
and 2 / = Qi, we guarantee, as anticipated, a correct low- 
energy description. The Canterbury group [HI ED] El] 
demonstrated that for Stieltjes and meromorphic bivari¬ 
ate functions, the convergence of C'^^'^(Qi, Qi) is guar¬ 
anteed for J = —1,0, properties exploited here [28]. 

The first element on the CA sequence reads 


cUQlQl) 


_Qo,0_ 

1 + ^i,o(Qi + Qi) + ^i.iQiQi 


(4) 


where Bose symmetry is already implemented (5^ = bj^ 
and 60.0 = 1 without loss of generality. Eq. Q repro¬ 
duces the high-energy behavior when one photon vir¬ 
tuality is set to zero, the well known Brodsky-Lepage 
limit E], property 5) above. Knowing the Taylor ex¬ 
pansion of the F'(Qi,Q 2 ), Eq. Q would be unique: 
00,0 = ^^( 0 , 0 ) is determined from the r( 7 r° —>• 77 ) 
through the relation ( 47 rQ:)^m^ai q = 647rr(7r° —>• 77 ); 
61^0 is the slope of the single virtual 7 r°-TEE and 61^1 
is related to the doubly-virtual slope. 

The second element on the CA sequence with the ap¬ 
propriate high-energy behavior results in 
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Cl{QlQl) = 


00,0 + oi,o(Qi + Q 2 ) + oi.iQiQi 


1 + h,oiQi + Ql) + KiQlQl + ^2.o(Qt + Qi) + &2.i(QtQ: 


Q\Qi)~+b2aQW2 


(5) 


and demands the knowledge of four coefficients belonging 
to the double virtual sector (ai.i, 62,1 and 62,2) as 

well as the curvature and third derivative of the single 
virtual TFF to match oi.o and &2,o- Obeying the con¬ 
vergence properties of the CA, C\{Q\,Q\) approximates 
better the TFF than Qi)- The difference among 

them can be taken as a way to estimate the systematic 
error on the approximation sequence [28) . 

Experimental data for the doubly virtual F{Q\,Q^) 
are not available yet and we cannot extract all those 
terms from them. To yield a result with our method 
we need to find a compromise. On the one hand, the 
OPE tells us that ^ Q~'^ and im¬ 

plies = 0 in Eq.Q and 62,2 = 0 in Eq.®. If the 
OPE is fulfilled, no experimental information about the 
doubly virtual TFF is required up to the second element 
on the sequence, the 02(Qi,!?!)- On the other hand, 
if we impose the factorization approach to Ci{Q\,Q 2 ), 
which means Ci{Q\,Q\) = x PiiQ^), we would 

find bi I = Q. In this second scenario, the high-energy 
tail of the TFF would behave as instead. How¬ 
ever, the yPT study performed in Ref. [21] favors the 
factorization approach at very-low energies for the TFF, 
indicating that the role of the high-energy tail may not 
be that relevant at the low scales which are crucial in 
our calculation m- In summary, in order to remain on 
the conservative side, we assume for the C'2’(Qi,Q2) 
proximant the —yet unknown— 61 1 value to fie within 
the reasonable 0 < bi^i < 6 ^ q range, which represents 
a compromise between the sought low-energy descrip¬ 
tion and the high-energy QCD constraint. The lower 
limit coincides with the OPE constraint and prevents as 
well the Ci{Q 1 ,Q 2 ) to have a divergency on the space¬ 
like region; the upper one is the value implied by fac¬ 
torization. In calculating the BR we do not use (|^, 
but (|^, including in addition the SM Z-boson contri¬ 
bution (of order —0.02 x 10“®). With such parame¬ 
ters, and taking the 6,r value obtained in |25j . we obtain 
BRsm( 7J‘° —>■ e“'"e“) = (6.20 — 6.35)(4) x 10“® where the 
two main numbers come from ranging 61.1, and the error 
comes from r(7r° —>■ 77) and b^ uncertainties. 

The next element of the sequence, the (Qi, Qi): 
demands (c.f Eq. (§) not only imposing the 

limQ2_>oo E((5^, Q^) ~ Q~^ but also its ex¬ 

act coefficient. For that reason, to extract 017 
and 627 in (j^, we take \imQ 2 ^^ F{Q'^,Q'^) = 

I§ 2 + 0 {Q-^)^ where = 0.20(2)GeV2 

from Ref. |31j . The remaining parameter in ([^, 
the 617, is associated to the low-energy expansion 
Tp7*7* = Fpj^{l + ^{Q\ + Q\) +... + + - -Oj 

i.e., with the yet unknown slope of doubly virtuallity 
«7r:i,i- To obtain 617, we determine a,r;i,i similarly 


as we did for the case before which is, 

after matching all the other parameters we let this one 
range freely with the only constraint of not developing 
divergencies in the space-like region. In this respect, 
we find 1.926^ < a,r;i,i ^ 2.076^. Experimental data 
will ultimately decide on this issue since knowing both 
the low-energy parameters and the OPE we will reach 
higher order in the Q 2 ) sequence and the 

discussion about OPE versus factorization will be 
unnecessary. Notice that the (Qf, Qi), even though 
has its parameters fixed by OPE constraints, still fulfills 
the factorization approach at low energies as discussed 
in Ref. |22|- 

With all these constraints, we determine a new SM 
range 

BRsM(7r° ^ e+e“) = (6.22 - 6.23)(4)(2) x 10“® , (6) 


where the two main numbers come from the ranging 
of 0^77. As before, the error ±4 • 10“® comes from 
F(7r° —>■ 77) and b^^ uncertainties, and the ±2x 10“® from 
the evaluation of the systematic error of our approxima¬ 
tion. Such error is obtained after comparing the differ¬ 
ence between the C'°((5i,Q2) ^od C 2 (Qi,Q 2 ) [13 HE]- 
The result in Eq. ® can be summarized as BRsm( 7>'° —>■ 
e“'"e“) = 6.23(^(3)^x 10“® with the systematic and the 
spread in Eq. ® combined linearly. This results reduces 
the error on the previous SM determination by more than 
50%, which demands now including the Z boson contri¬ 
bution and makes no approximation on the calculation of 
the loop integral. Given that our final error is dominated 
by the input errors and not by the observed convergence 
of our approximants, we anticipate that the next element, 
the Q 2 ), would not provide any valuable informa¬ 

tion. In this respect, we believe that the expected exper¬ 
imental measurement of the slope of doubly virtuality 
will not make a dramatic change with respect to Eq. 
beyond reducing the range quoted (unless we also find a 
surprise there). 

For completeness, let us indicate that if we would have 
considered approximation (§ instead of the full result 
derived from ([^, we would have obtained (6.17 — 6.18) x 
10“®, while omitting the Z boson contribution would lead 
to (6.24 — 6.25) X 10“®. Eq. Q still represents a devi¬ 
ation of the measured BR of about 3.2tT (or l.Scr with 
RG in [171 118)1. which we think is unsurmountable in a 
model-independent way with our present knowledge of 
the Standard Model. Therefore, to eventually improve 
on the situation, experimental data would be, as said, 
required. 
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IV. IMPACT OF THE KTEV MEASUREMENT 

ON THE HLBL CONTRIBUTION TO {g - 2)e,^ 

CAs are flexible enough to be able to reproduce, for 
the first time, the KTeV measurement —which translates 
into tuning 01^,611 in Eq.([^ — if our previous scenarios 
(i.e., factorization and high-energy QCD constraints) are 
momentarily put aside, a procedure which does not spoil 
any property from Fade Theory m- 

This result is shown in Fig. [^as a purple band when 
including the new RC, for which —396^ < a,r;i.i < 

For completeness, we show as well what would have been 
achieved without the new RC as an orange band (which 
corresponds to —3706^ < a 7 r;i,i < —1006^). While the 
very low-energy region is not very different from our re¬ 
sult Q (blue band), KTeV seems to imply a strongly de¬ 
creasing TFF at large Q^. This effect could be achieved 
with a very slow converging OFF as indicated by the 
large (5^ > 10 GeV^ value we obtained through match¬ 
ing aiq and bi^i to KTeV. This feature has never been 
observed before. 

A first glance on Fig. [^reveals that, even without high 
precision data at a given (30% or even 50% exper¬ 
imental errors), these 3 scenarios could be easily dis¬ 
tinguished, while improving our result ([^ would require 
higher precision, around (10 — 20)%. Such data may be 
fitted through CAs, which being able to accommodate 
the high-energy constraints from QCD as well, would al¬ 
low to reconstruct the TFF from = 0 to 00 in the 
space-like region. 

Given that nowadays the KTeV measurement repre¬ 
sents then the only (indirect) source of experimental in¬ 
formation on bi I, as an amusement, we can use the re¬ 
sults from the previous paragraph and explore the im¬ 
pact of KTeV on the hadronic Light-by-Light contri¬ 
bution to {g — There, the 7r° yields the lead¬ 
ing contribution, and thereby the 7r°-TFF 

plays a central role in the calculations [32] • Taking the 
range given by the purple band, we find for the (off- 
shell) 7r°-exchange = 2.0(3) x 10“^"^, where 

the error accounts for the band’s width, (to compare 
with 3.0(3) X 10“^^ |321)- For the (on-shell) 7r°-pole 
we obtain = 1.8(2) x 10“^^ (to compare with 

2.6 X10-14 jsni). 

Given that = -105(81) x IQ-^^ (32(33], this 

effect is negligible though in the correct direction. This 
contrasts with its muon counterpart, which in addition is 
more sensitive to hadronic physics. There, the deviation 
a^P — = -1-290(90) x lO-n [32 ES] has opposite sign. 

For the (off-shell) 7r°-exchange we obtain = 

40(9) X 10-11 (to compare with 72(12) x lO-n (32). For 
the (on-shell) 7r°-pole we obtain = 36(7) x lO-n 

(to compare with 58(10) x lO-n [53]). These shifts rep¬ 
resents twice the foreseen experimental accuracy (16 x 
10-11) in the future {g — 2)^ experiments projected in 
Brookhaven and J-FARC (33] [3S] indicating that the 
current precision of the SM error on the {g — 2)^ would 


be underestimated if the > e^'e- would be taken 
into account. Taking KTeV results without the latest 
RG would accentuate the differences indicated above. A 
thorough application of our approach to (g — 2) will be 
discussed elsewhere. 


V. NEW PHYSICS CONTRIBUTIONS TO 


As anticipated in the introduction, tree level contribu¬ 
tions from new physics may be relevant for this process 
as well. In our case, only pseudoscalar (P) and axial 
(A) channels are relevant, since interactions such as lep- 
toquarks may be Fierz-rearranged [35], then, only pseu¬ 
doscalar and axial channels will appear. The following 
Lagrangian 


^ (/^T's/) + 2mfcf P, (7) 


where g,rnw are the standard electroweak parameters, 
and Cj are dimensionless couplings to the fermion / = 
{u,d,e}, yields an additional term modifying Fq. ([^ 


A{q^) ^ A{q^) + 


V2GfF^ 




{X^ + Xn, ( 8 ) 


where Gp is the Fermi coupling constant, the 
pion decay constant, = cf — c^) and X^ = 
cf (cu ~ c^) /(I ~ Wp/m^o)- As an illustration, the 
boson = —c^ = — 1, = 0) shifts A{q^) by —0.3%. 

Similarly, the “dark” Z model in [16], would reduce the 
amplitude by (—0.3%)^^, with ^ 1 [37]. Therefore, 
electroweak-like effects are unlikely to be important and 
more general NF approaches would be required when 
considering axial interactions [14]. Besides, pseudoscalar 
contributions appear in extended Higgs sectors, such as 
the supersymmetric model in Ref. m- Their effects in 
our decay may become large enough if mp ~ 171^^0 as 
may be inferred from the pole in X^, though existing 
constraints on nip discard this scenario. Since many 
phenomenological constraints to these simple NF sce¬ 
narios exist [37] (dark photon searches, {g — 2)ep, elec¬ 
troweak physics), a suitable scenario represents a chal¬ 
lenging study which is beyond the scope of this article. 


VI. CONCLUSIONS 

In this work, we presented a model-independent ap¬ 
proach to describe symmetric bivariate functions based 
on Canterbury approximants. This method allows a 
model-independent study of the doubly virtual 7r°-TFF 
in the space-like region which incorporates the low- and 
high-energy QCD constraints. These are essential prereq¬ 
uisites for calculations such as 7r° —)■ e+e- or . 

We have predicted BRsM(7r° —>■ e+e-) = 6.23(4)(3) x 




6 



CN 


9 






[GeV^] 


[GeV^] 


FIG. 3: Left panel: normalized TFF assuming Q\ = Q 2 = ■ Right panel: normalized TFF assuming Q\ = 0.5 GeV^. 

Black solid line indicates the factorized TFF. Upper (blue) band shows our CHQ^jQV) estimation with 
1.925^ < a.n-;i 4 < 2.076^. Lower (orange) band reproduces the KTeV measurement within Ict. Middle (purple) band considers 

KTeV measurement with the new RC 


10“®, still 2cr off the experimental result and found that, 
unless New Physics are present, the KTeV result on 
TT^ —)■ e“'"e“ implies an unexpected behavior for the dou¬ 
bly virtual TFF. This would produce a large shift for 
, of the order of projected experiments for mea¬ 
suring {g — 2)^. The effect of CAs in 7r° —>• 4£ decays 
is studied in [3S]. Possible effects may appear as well in 
ep(en) elastic scattering, which would be relevant for the 
proton radius puzzle. 

We conclude that the current situation demands an 


experimental measurement of (Qf, well as 

a new —>■ e“'"e“ determination. 
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